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Abstract: Baxter's Q-operator for the quantum transfer matrix of the XXZ spin-chain 
■ is constructed employing the representation theory of quantum groups. The spectrum 

of this Q-operator is discussed and novel functional relations which describe the finite 
. temperature regime of the XXZ spin-chain are derived. For non-vanishing magnetic field 

the previously known Bethe ansatz equations can be replaced by a system of quadratic 
equations which is an important advantage for numerical studies. For vanishing magnetic 
field and rational coupling values it is argued that the quantum transfer matrix exhibits a 
\^ ' loop algebra symmetry closely related to the one of the classical six- vertex transfer matrix 

at roots of unity. 

^— > • 

d • 1. Introduction 



In this work we present new identities for the description of the spectrum of the quantum 
transfer matrix for the XXZ spin-chain with non- vanishing external magnetic field h, 
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tfxxz = ^ E K^+i + + A(*Mf 1 - !)} - \ E °* ■ (i- 1 ) 

Here A = (q + q~ 1 )/2 is an anisotropy parameter which is assumed to be real. This spin- 
chain serves as a prototype model for other more complicated integrable systems. It is an 
important toy model for the exact computation of physical quantities such as magnetic, 
electric or thermal conductivities. In this context the study of the finite temperature 
behaviour of the spin-chain is of crucial importance. One method to achieve this is the 
so-called quantum transfer matrix (see e.g. O), 

t(z;w) = Trq a R m {zw)R\f_ 1)Q (w / z) ■ ■ ■ R 02 (zw)R\f (w/ z), q a := e^ 2 (1.2) 

in terms of which the partition function of the XXZ spin-chain can be expressed 

Z L = Tr e~ 0HxKX = lim Tr t(z = 1; w = e' 13 ' /N ) L , = 8(q — q^ 1 ) . (1.3) 
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Let us explain the various objects appearing in the definition. The variable (5 > denotes 
the inverse temperature of the system and for convenience we have introduced a "twist 



angle" a = /3h/2\nq with h being the magnetic field in (1.1). The quantum transfer 
matrix is built out of the well-known six-vertex .R-matrix which acts on the tensor product 

C 2 ®C 2 , 



R=\ ®cr + cr + <g> a cr + + c a + <g> a" + c' a' <g> a + . (1.4) 

\ b I V a I 

Here the parametrization of the Boltzmann weights a, b, c, c' is chosen as follows 1 

! , (1 ~ g)g 1 / n fr\ 

a = 1, o=- 5-, c=- 5 , c = cz . (1.5) 

1 — zg z 1 — zg z 



The upper index i <g> 1 in ( |1.2| ) stands for transposition in the first factor. We recall that the 
six-vertex i?-matrix gives rise to a classical statistical mechanics system which is described 
by the classical six- vertex transfer matrix (as opposed to quantum), 

t 6v (z) = Tr R 0L (z) ■ ■ ■ R 01 (z) . (1.6) 

This classical physical system is connected with the above quantum spin-chain through the 
relation (we set temporarily h = 0) 



#xxz = (<?-<? 1 ) z-^-lnt ev (z) 



(1.7) 

=i 



The prefactor in the last equation explains the introduction of the renormalised temperature 
variable in (|ll|). As an immediate consequence of ( |1.7| ) one has the following identity 
for the density matrix 

lim ft 6v (l)- 1 t 6v (e^'/iV ) N ) JV = e -/3H xxz (L8) 
This rewriting is inspired by the Trotter formula in the path-integral formalism of quantum 



field theory and N is referred to as Trotter number Q. The last expression (1.8) can be 



conveniently expressed in terms of the quantum transfer matrix ( |1.2[ ); see the review |2[ 
and references therein for details. Note in particular that for this construction to work 
the Trotter number N has to be even and we shall work throughout this paper with the 
convention 

N = 2n . (1.9) 

In the thermodynamic limit when the physical system size L tends to infinity the 
partition function (|l,3|) can be approximated by the largest eigenvalue A^r of the quantum 
transfer matrix (|L 



L > 1 : Z L = lim Tr r(l; e' 13 '^) 1 rs lim (Aiv) L . (1.10) 



1 Below we will compare the results in this paper with the known properties about the quantum transfer 
matrix as they can be found in e.g. To this end it is helpful to identify the definition of the Boltzmann 
weights in on page 11, equation (2) with ours by setting z — e 11 ™ and q = e 17 . Note that on page 16 in 
pj] a rotation in the complex plane is performed replacing w — ► iv; see equations (23) and (24) therein. 
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Thus, the introduction of the quantum transfer matrix allows one to efficiently in- 
vestigate the finite temperature regime. Instead of having to compute multiple excited 
states and energies of the quantum spin-chain ( |1 , 1|) at finite temperature, one only needs 
to compute a single eigenvalue and eigenstate of the quantum transfer matrix (1.2). Similar 



simplifications occur also in the computation of finite temperature correlation functions; 
see for instance ||. 

There are some technical complications, however. We refer the reader for the following 
statements to Q and references therein. The algebraic properties of the quantum transfer 
matrix resemble closely the ones of the classical six-vertex transfer matrix ( |1.6|) and as a 
result one can compute the eigenstates and eigenvalues of the quantum transfer matrix by 
similar methods as in the classical case, i.e. the algebraic Bethe ansatz. Via this route one 
obtains expressions for the eigenvalues of (|l]^). In particular the largest eigenvalue A^r can 
be cast into the form 

. e^(/){u-i)Q(u + 2i) +e~ S T<f)(u + i)Q(u-2i) 

An{u) = i\r , (l-U) 

Q(u)[sinh ^(u — 2i + it) sinh ^(u + 2i — it)] 2 

where we have set z = e 7 ", w = e _/3 / N = e _47r , q = e* 7 and introduced the functions 

n 

(f)[u) = [sinh %(u-i + ir) sinh ^(u - 2i + ir)] n , Q(u) = U sinh %(u - Uj) . (1.12) 

j= i 

The quantities Uj are solutions of the following system of non-linear equations 

p±j) ^a^) „ (113) 

4>[Uj — 1) Q{Uj — 2i) 

Compare with equations (25-30) in Q. The analytic solutions to the last set of equations, 
known as Bethe ansatz equations, are not known. Furthermore, in the Trotter limit N — > 00 
the distribution of the Bethe roots Uj remains discrete and cannot be approximated by 
continuous density functions as it is the case for the classical transfer matrix ( |1.6| ). Instead 
one has to rely on the numerical solution of a non-linear integral equation. The derivation 
of this integral equation as well as other properties of the quantum transfer matrix are 
based on numerical studies of the above Bethe ansatz equations ( |1.13| ) . 

With this in mind it is worthwhile to study alternative methods of deriving the spec- 
trum, and in particular the largest eigenvalue, which lead to systems of equations simpler 
than the Bethe ansatz equations ( [1.13 ). The purpose of this article is to diagonalise the 



quantum transfer matrix employing Baxter's idea of an auxiliary matrix known as Q- 
operator [||. The eigenvalues of the Q-operator give the Q-function in ( [1,12 ) and we will 



derive a set of functional relations analogous to the ones first obtained in the context of 
conformal field theory ||. The case of the finite XXZ spin-chain has been investigated in 
and |8|, 0, In particular, we will follow closely the treatment for the twisted XXZ 



spin-chain given in [1C] 



In order to construct the Q-operator we will use the representation theory of quantum 
groups. The main results for h 7^ in ( |1.1| ) are novel identities for the spectrum of 
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the quantum transfer matrix and a simpler set of equations which imply the Bethe ansatz 
equations ( [1.13 ) but are quadratic instead of order N. One of the main results in this paper 



is that the largest eigenvalue A^r of (1.2) can be expressed in terms of one polynomial 



Q+(z) = £e+(-z)* 



(1.14) 



k=0 



whose coefficients (with = 1) solve the following system of quadratic equations 



E 

k+l=m 



n \ I n 



k V I 



(wq) 



k-l 



E 

k+l=m 



sinh[^ - i 7 (Jfe - I)} 



sinh[/3/i/2] 



(1.15) 



Here the summation convention in ( |1.15 ) is to be understood as follows. First fix the 
variable m in the allowed range m = 1, ...,N — 1 and then sum in ( |1.15| ) over all possible 
values for k, I such that k + I = m. Thus, one obtains in total N — 1 coupled quadratic 
equations. For real z the largest eigenvalue of the quantum transfer matrix is then given 
by 



A N (z;w) 



Q + (zq~ 2 )Q-(zq 2 ) - e-P h Q + {zq 2 )Q-{zq- 2 ) 



smh\3h/2][(zwq - l)(z/wq - l)] N l 2 
where Q~ is the reciprocal polynomial of Q + , i.e. 

n + 

Q~{z) = E ^HO* = ^ n Q + (^ 1 )/e+ • 



(1.16) 



k=0 



C-n 



In the critical regime, g £ S 1 , the coefficients et obey the additional constraint 



(4r 



(1.17) 



There are in general many solutions to the equations ( 1.15| ) (similar as there are multiple 



solutions to the Bethe ansatz equations) describing a subset of the spectrum of (1.2). The 
largest eigenvalue of the quantum transfer matrix Ajv appears to be always among them; 
this has been numerically verified for N = 2, 4, 6, 8, 10, 12, 14, 16 2 and for all of these cases 
the total number of solutions to (|1.15 ) is found to be 2 n . 

The identities ( |1.15 ) and ( |1.16 ) are special cases of more generally valid expressions 
for the spectrum of the quantum transfer matrix; see equations ( [4.2; ) and (|4.3| ) in the text. 
We discuss these identities in the context of higher spin quantum transfer matrices and the 
associated fusion hierarchy. This is motivated by the definition of a trace functional which 
might be relevant for the computation of correlation functions. Underlying the derivation 
of the main result is a concrete operator construction for which is important to extract 
further information on the spectrum and also applies to the case of vanishing magnetic 
field, h = in flnj), albeit one has then to restrict the deformation parameter q to a root 



2 It appears that there is a difference between n = N/2 odd and even; compare with the footnote in j3|, 
Section 2.6. For n odd the largest eigenvalue might have less than n Bethe roots. Nevertheless, in the cases 



N = 2, 6, 10 we have numercially verified that (1.15) and (1.16) still hold true 
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of unity, q l = 1, I > 2. In this case, we argue that the quantum transfer matrix exhibits 
a loop algebra symmetry just as the classical six-vertex transfer matrix albeit in a 
different representation. 

The outline of the article is as follows. In Section 2 we connect the structure of 
the quantum transfer matrix to the representation theory of the quantum affine algebra 
U q {sl2)- In particular, we discuss a specific dual representation and the corresponding 
L-operator. This will yield the commutation relations between the quantum monodromy 
matrix elements corresponding to ( |1.2[ ) and the Chevalley-Serre generators of U q {sl2) which 
allows us to prove the loop algebra symmetry U{sl2) of the quantum transfer matrix at 
roots of unity and for vanishing magnetic field. As a preparatory step for the discussion of 
the Q-operator we introduce the quantum fusion hierarchy, i.e. the higher spin analogues 
of the quantum transfer matrix. 

In Section 3 we construct the Q-operator for the quantum transfer matrix and discuss 
its properties and functional relations. The proofs can be found in appendices A and B. 
For the construction of the Q-operator one has carefully to distinguish between the case of 
generic q and q a root of unity. In the former case the auxiliary space of the Q-operator 
is infinite-dimensional and one needs to introduce a boundary parameter (the external 
magnetic field h) in order to ensure convergence. At a root of unity the auxiliary space is 
finite-dimensional and the construction of Q then also applies to the case h = 0. However, 
some of the functional relations which hold true at h ^ then cease to be valid. 

Section 4 is devoted to a special Q-operator functional equation, the Wronskian rela- 
tion, which suffices to compute the spectrum of the quantum transfer matrix and implies 
the Bethe ansatz equations. At the end we discuss some special solutions which contain 
the largest eigenvalue and are based on numerical evidence. 

Section 5 contains the conclusions. 

2. The Quantum transfer matrix and representation theory 

As a preparatory step for the construction of the Q-operator for the quantum transfer 
matrix let us first analyse the definition ( |1.2[ ) from a representation theoretic point of view. 
This will enable us to define the corresponding L-operator from which the Q-operator will 
be built. 

Recall that the six-vertex i?-matrix is an intertwiner of the quantum affine algebra 
Uq^sfa) with respect to the tensor product of the two-dimensional evaluation representation. 
The quantum affine algebra U q {sl2) is generated from the Chevalley-Serre elements subject 
to the relations 

qhtqhi =q hiqh, q h ^e jq ~ h ' = q A '^ ej, q h ' f^ = q~ A ^ f v fe, /,•] 
and 

X Xj XjX<i | ^3 j q X i X j X ^ XjXj^ — 0} X — €■ ^ J~ 



q - q- 1 

(2.1) 

(2.2) 
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Here the indices i,j take the values 0,1 and Aij is the Cartan matrix of 5/2- The evaluation 
homomorphism ev z : U q {sl2) — > U q {sl2) defined by 

ev z (e )=zf, ev(f )=z- 1 e, ev(q h °) = q~ h (2.3) 

and 

ev( ei )=e, ev(h) = f, ev(q h ') = q h . (2.4) 

An evaluation representation is now obtained by combining the evaluation homomorphism 
with any finite-dimensional representation of U q (sl2), in particular we can choose the two- 
dimensional, spin 1/2 representation in terms of Pauli matrices, 

7r(e)=a+, tt(/) = <7- ir(q h ) = tf* . (2.5) 

The six-vertex R-matrix then intertwines the tensor product representation ir z (g) ix\ with 

TT Z = 7T O eV z . 

Given any representation p : U q {sl2) — * EndV over some finite-dimensional vector 
space V we define the following representation over its dual space V*, 

p* : U q {sh) — > EndF*, (p*(x)v*,w) := (v *, p^ 1 (x))w) , x € U q (sh), w 6 V . 

(2.6) 

Here 7 is the antipode which is defined on the Chevalley-Serre generators as follows, 

7 (ei) = -g^e*, 7 (/i) = -/i<?\ 7^) = f 1 * (2-7) 
7-Hei) = -e iq ~ h % 7 - 1 (/ f ) = -^/i, 7 _1 (^) = ^ • (2-8) 

If we canonically identify with V the representation in terms of matrices is simply 
given by 

p*(x) = (p( 1 - 1 (x))f, xeU q (sl 2 ) (2.9) 

with t denoting the transpose. Note that the definition ( |2.7| ) is compatible with the fol- 
lowing choice for the coproduct 

A(ei) = ei ® 1 + ® e i5 A(/i) = /< ® g"' 1 * + 1 <g> A(g^) = g^ ® q h > . (2.10) 

Setting p = it, the two-dimensional evaluation module, we are interested in finding the 
intertwiner 

L* A*(x) = A* op {x)L* with A* = (1 ® vr*)A and L * G ^(^2) ® EndF* . (2.11) 

Explicitly, the coproduct relations for the Chevalley-Serre generators read 

A*(ei) = ei®l-q 1+hl ®o-~, A* p (ei) = e x O g~ CT * - g 1 ® <r~ , (2.12) 

A , (/i) = /i®^-1~ 1 l^ + , A;(/ 1 ) = / 1 ®1-<|- 1 -' !i ® ( 7 + , (2.13) 

and 

7r*(ei) = -(^ + g- CT2 ) t = -g ^, tt*(/i) = -(^V-) 4 = -g~V+, (2.14) 

7r*(e ) = -(a- q °y = -q a+, 7r*(/ ) = -(g~ CT V+)* = -^V" . (2.15) 
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We will see below that the intertwiner (|2,11 ) specializes to the i?-matrix used in the defini- 
tion of the quantum transfer matrix ( |1.2[ ). Solving the above intertwining condition ( 2.11| ) 
for an evaluation module with central charge zero, i.e. q h ° = q~ hl , we find 

/ ft l+l h 1+ l \ 

t*( \ i zq 2 ~Q 2 W -Q)eiq 2 \ 
L(z)=[ hl+1 • (2.16) 

\z(q - q)q 2 f 1 zq 2 - q 2 J 

For comparison and in order to keep this article self-contained recall that the conventional 
L-operator reads 

' ~"~ Ll Z( "^ )q U 1 ) (2-17) 
\{q-q Jeig 2 Z q 2 - q 2 J 

In light of the following identity for the universal R-matrix (1 ® j~ 1 )R = R 1 it would be 
more natural to use — z _1 L*(z) as intertwiner for the dual representation ir*. However, we 
wish for later purposes to keep the L*-operator polynomial in z instead of z . Evaluating 
this intertwiner in the two-dimensional spin 1 /2 representation (|2.5|) yields 

= (ir®l)L*(z) = ^_ia«t = [i?21 ^-i)]i®t (2.18) 
zq 1 — q 



which is one of the R- matrices used in the definition of the quantum transfer matrix (|b^) . 
Thus, we conclude that one lattice row associated with the monodromy matrix of the 
quantum transfer matrix corresponds to the following quantum group module 



71% 



gjtW, gn<M = 7T* ®7r tt -i---® 7r;®7r m -i . (2.19) 



N 



From this we immediately deduce that the quantum transfer matrix t(z; w) block decom- 
poses with respect to the following alternating spin operator, 



1 N 

t(z;w),S A } = 0, S A :=-J2H k °k, (2-20) 



fc=l 

since the quantum monodromy matrix 

r.../~\ d c„...M3*®l/ 



R 0N (zw)R^_ 1)0 (w/z) ■ ■ ■ R 02 {zw)R\^(w/z) = I ^ 



(2.21) 



is an intertwiner with respect to the tensor product ( 2.19[ ). More generally, we have 

[A, q Hl ] = [D, q Hl ] = 0, q Hl Bq~ Hl = q~ 2 B, q Hl Cq~ Hl = q 2 C, q H ' = q 2SA . (2.22) 



Denoting the Chevalley-Serre generators acting in quantum space ( 2.19 ) by capital let- 
ters, {E\, Eq, F\, Fq, Hq = —Hi}, one finds the analogous commutation relations as in the 
classical case; see equations (14-15) in [p^ . The difference between the relations for the 
classical and the quantum transfer matrix is purely in the explicit form of the quantum 
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group generators which is fixed through the identification of the quantum group module 
aJt^in pl9|) . For instance, employing ( |2~10|) and (^), (glj) one has 



N 



Ei = J>g * II 4 k , ^ = (-1)" (2.23) 

fc=i v<* / 

*i=£ etg^o-r* f n 9"* JflJ ) > (2-24) 
fc=l V >fc / 



and so forth. Following the same line of argument as presented in [|T3|] (c.f. equations 
(23-24) therein) one then proves that for zero magnetic field h = and q being a primitive 
root of unity of order I the quantum transfer matrix Ql.2j ) enjoys a loop algebra symmetry 
U(sl2) in the sectors 

25yi = 0mod^. (2.25) 

The Chevalley-Serre generators of the loop algebra U^sfa) are obtained from the restricted 
quantum group (compare with the discussion in J14| |) via taking the following limit from 
generic q' to the root of unity value q, 

q'^q q — q I «/2, if £ is even 

Analogous expressions hold for the remaining generators. All of the Chevalley-Serre gen- 
erators, {E^f\ Eq , F± ,Fq}, commute in the commensurate sectors ( 2.25| ) with the 
quantum transfer matrix, e.g. 

[r(z;w),Ep] = for h = 0, q £ = 1 and 2S A = Omodf . (2.26) 

This result for the quantum transfer matrix is analogous to the loop symmetry of the 
classical transfer matrix first discovered in [12|, albeit via a different proof. 

2.1 Transformation under spin reversal 

In order to show the existence of two independent solutions to the TQ-equation, we discuss 
the behaviour of the quantum transfer matrix under spin-reversal. From the elementary 
identities 

(1 ® a x )R(z)(l ® a x ) = (a x ® z~^)R(z){a x ® z^) 
(1 ® a x )R*{z){l ® tr") = (cx* ® z^ir^o* ® 

we infer that the quantum transfer matrix (|1.2| ) transforms as 

r a (z; io) = r_ a (z; to) (2.27) 

under the involution w 9t with 



TV * * 



9\ 



\[a x and S* = ^^. (2.28) 



i-i ■ 2 

i=i 
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In addition, the quantum transfer matrix obeys another identity. First we observe that the 
following equations for the Boltzmann weights ( |1.5| ) hold true, 



a z -i = 1, 6,-1 = l/b zq -2, Cg-i = -q c' zq - 2 /b zq - 2 , c£_i = -q 1 c zq -2/b zq - 2 . (2.29) 
Employing the identities 

R(z- r ) = — !— (a x ® {-q)-^)R{zq- 2 ) lm (o- x ® (-q)^) 
R*{z' X ) = {a x ®{-q)^r)R*{zq 2 ) mt {u x ®{-q)-^) 

z -l q -2 

we then easily find the expression for the transpose of the quantum transfer matrix, 

1 „.,-l\ _ TU. „—a ° z ®^Jl -,-!„.. D* D_ . 1„,.— 1 

-2\t 



_ r a (z, Wg - , 

where we have used that [SA,T a (z,w)] = 0. 
2.2 The quantum fusion hierarchy 

The quantum transfer matrix (|1.2[ ) commutes with an infinite family of higher spin transfer 
matrices. In close analogy with the classical six-vertex model we define for d £ N the 
following family of transfer matrices, 3 

T (d-i)( z . w)= Tt q ahm L N (zw)L* N _ l (z/w)---L 2 (zw)Ll(z/w), q a = e^ , (2.31) 

where the representation in auxiliary space has been replaced by the spin (d — l)/2 evalu- 
ation module 7r( rf_1 ), 

^ d ~ X \e) \k) = [d-k] q \k-l), TT^Cf) \k) = [k+l] q \k+l), ^ d - l \q h ) \k) = q d - 2k - 1 \k) , 

(2.32) 

where the index k labelling the basis vectors of the representation takes values in the set 
k = 0, 1, 2, d — 1 and we set e |0) = / \d — 1) = 0. Our motivation for introducing ( 2.31| ) 
is twofold. They are natural objects to consider from a representation theoretic point of 
view and we will encounter them when deriving functional relations for the Q-operator in 
the subsequent section. The other reason is their extension to complex dimension d E C 
which is closely related to the trace functional used in recent formulations for correlation 
functions ||TTf| . 

Setting d = 2 we identify -n^ 1 ) = it in (2.5) and recover the quantum transfer matrix 
via the relation 

= (-z/w)-»tW(z;w) = tW(z;w) 
' (zwq — q 1 ) n (wq/z — q x ) n (zwq — q 1 ) n (z/wq — q) n 



3 In the following we shall often suppress the explicit dependence on the temperature parameter w. 
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The spin representation yields the quantum determinant, 

t(°\z;w) = (zwq^ — q~^) n (zq~^ /w — q^) n 

= (zw - q- l ) n {z/w - q) n = {zwq - l) n {z/wq - l) n (2.34) 

Similar to the classical six-vertex model, the higher-spin quantum transfer matrices r" 
satisfy a functional relation known as the fusion hierarchy, 



^ d - 1 \zq d )T^\z)=T^\zq- 1 )^ d - 2 \zq d+1 )+T^\zq)T^{zq d - 1 ) (2.35) 
rollan 

the exact sequence 



which is a corollary of the decomposition of the tensor product 7r^ d ^ ® nz described by 

zq 



- tt - ^ S vr^L - . (2.36) 

Since the auxiliary spaces in the quantum transfer matrices are the same as in the classical 
case, we can use the same representation theoretic results to derive all relevant functional 
relations. What changes in the transition from "classical" to "quantum" are the coefficient 
functions which appear in the respective functional equation. For instance, the coefficients 
T^°\zq~^), r^°\zq) follow from the identities 

(vr^- 1 ) ® l)L 13 (zq d )R 23 (z)(l ® 1) = (z - l)(t ® l)^ 2 ) ® l)L(zq d+1 ) 

(p (8 1)(tt^ ® ljLia^i^^) = [zq - q' 1 ) (tt^ ® ® 1) (2.37) 



and 



(p® ^(vr^- 1 ) ® 1)L* 3 (V)^23(^ = («- 1 -l)(7rW®l)L*(2g d - 1 )(p®l) (2.38) 



where the maps 



and 



i : •= — >- [<i — /c — 1], |fc) <8> |1) - g^*" 1 ^ + 1] 9 I*: + 1) ® |0) (2.39) 



f : JJ^k] l*> ® l°> "» ( 2 - 4 °) 



are the inclusion and projection map in the exact sequence ( |2.36| ). 
3. The Quantum Q-operator 



After outlining the derivation of fl2.35 ) we now turn to the Q-operator and apply the same 



strategy as in the case of the fusion hierarchy. The main difference lies in the fact that for 
the definition of the Q-operator we need to introduce an m/miie-dimensional evaluation 
module when q is generic (i.e. not a root of unity) Q ||, 



p+(e ) \k) = Z \k + 1) , P + (q^) \k) = p + (q-^) \k) = r V*~ 1/2 
s + 1 — q 2k — sq~ 2k 

(q-q- 1 ) 



P + {ei)\k) = ,„ * \k-l), p+(e 1 )|0> = 0, r,s,zeC. (3.1) 
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Here s,r are free parameters characterizing the representation and z denotes the spectral 
variable as before. In connection with spin-reversal one also encounters the module 

p--=p + ouj with {ei,e ,q~,q 2 } -» {e ,ei,q 2 ,q~ } . (3.2) 

Note that in the limit s — ► we recover the g-oscillator representations used in ||. In the 
case that q is a primitive root of unity of order I we truncate the evaluation module p + by 
imposing the condition (compare with 0) 

= „ +M K-l>=0, ?=(/■ , i ;; iS ° dd . (3.3) 

I 1/2, it I is even 

The intertwiner corresponding to the quantum group module p + ® ir has been computed 
previously [0 || and reads, 

m=( Z " qh ^- q ~ h S ^-O^^eU^^EndV. (3.4) 

V (q — q ) e iq 2 zr q 2 ~ q 2 J 

In order to define a Q-operator for the quantum transfer matrix we now need to compute 
the intertwiner corresponding to the module p + ®tt*. The result is 

( zr q 2 - q 2 {q - q)eiq 2 \ ^ TT n v „ v , T ,* r „ ^ 

£ = _ hj+i hj-i € U q {sl 2 ) ® EndV . (3.5) 

\(9 -?)? 2 eo 2 - 9 2 / 

The last expression is derived from the coproduct relations 

A*(ei) = ei ® 1 - q 1+hl ® a", A* p (ei) = ei ® g _<7 * - g 1 ® a - , 
A*(eo)=e ®l-g 1+ ' io ® ( J + , A* p (e ) = e ® ^ - q 1 ® ^ + ■ (3.6) 

We define for r = 1 in p + the operator 

Q(z-s) = Trq^^QNizw^^iz/w) ■ ■ ■ £ 2 (zw)Z\{z/w), q a = e^ 2 , (3.7) 



where we have normalised (|3.5| ) such that Q is polynomial in the spectral parameter z. 
The specialization to r = 1 can be imposed without any loss loss of generality due to the 
identity 

Q(z;r,s) =r a - SA Q(z;r = l,s) . (3.8) 

Having stated the explicit definition of the Q-operator we now turn to its properties and 
the functional relations it satisfies. Since the auxiliary spaces of the quantum transfer 
matrix and the Q-operator are the same as in the conventional, classical six-vertex model 
the proofs for the statements made below follow closely the line of argument presented 
previously in || ||, [l0| . I shall therefore omit the proofs from the main text and refer the 
reader to the appendix for further details. 
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3.1 Factorisation 

One important result is that the Q-operator factorises into simpler operators as follows (we 
still assume a/0), 

Q(z; s) = Q(0; s)Q+ {z)Q~ (zs), (3.9) 



where 



Q+(z) = lim Q(0; s)~ l Q{z; s) (3.10) 



and the normalisation factor Q(0; s) at z = is easily computed to be 

1 



lim Q(z;s) = Tt q^ a ~ SA ^ hl = < 

2^0 p+ 



qot-S A — qSA-a ' 



g generic 



l_^-a) 

= ±1 



(3.11) 



qot—SA — qSA—a ' 



The identity (3.11) for generic g has to be understood as analytic continuation from the 
region of convergence. Note that the point a = 0, i.e. vanishing magnetic field h = 0, 
remains singular. This can be understood from the construction of the Q-operator since for 
generic q the auxiliary space given by the the representation ( p.l| ) is infinite-dimensional 
and twisted boundary conditions with an appropriate choice of a are needed to ensure 
that the trace in (3.7) is well-defined; compare with the discussion in H. The remaining 
operator Q~ is the counterpart of Q + and can be obtained by employing the following 
transformation, 

Q{z~ 1 ,w~ 1 ;s) = z~ n s% +Sa (w/q)- SZ 9\ Q{z/s,wq- 2 ; sf V\{w/q) SZ . (3.12) 

Here the operators S z , $K have been introduced earlier in ( p. 28 ) and use has been made of 
the identities 

2(z- 1 ) = - z -\(l®(-j-)-°\*)Z{ zq -^/ s Y®\l®a%-f q y z )s^ 
^{z- 1 ) = -z- 1 q- 1 {l®{-^y Z o x )2?{zq 2 /s) mt {l®o x {-^)- CTZ )s 1 -f - . 

The identity ( |3.12| ) then implies that up to an unimportant normalization constant we 
have, 

Q~(z,w) oc zf +SA (wq) sZ 9\ Q + {z~ l , W - l q- 2 ) t 9\(wq)- sZ . (3.13) 

Note that the action of the spin-reversal operator relates Q~ to the representation (|3.2j). 
Both operators Q ± have polynomial eigenvalues w.r.t. the spectral variable z. I shall 
denote these eigenvalues by the same symbol as the operators, 

n± n± 

Q ± (z) = l[(l-x±z)=J24(- z ) k > n ± = n^S A . (3.14) 

k=l k=0 

As we will see below the polynomial roots are two sets of Bethe roots. They can be 
(numerically) computed by employing a number of functional relations which are satisfied 
by the Q-operator. 
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3.2 Functional relations 

The best known functional relation is the generalization of Baxter's TQ-equation for the 
six-vertex model. This equation is obtained in the present construction for the Q-operator 
(which differs from Baxter's approach) by first deriving the functional relation 

Q(z; s)t^(z) = q a - SA T^(zq)Q(zq- 2 ; sq 2 ) + q SA ~ a r^ {zq- l )Q{zq 2 ; sq~ 2 ) (3.15) 

which is a direct consequence of the following decomposition of the tensor product of 
representations, p + (z; r, s) g> ir z , described by the exact sequence j?]] || 

— * p + (zq 2 ; rq" 1 , sq~ 2 ) <— > p + (z; r, s) (g) ir z — > p + (zq~ 2 ,rq, sq 2 ) —> . (3.16) 

Here r(°) is the quantum determinant introduced in fl2.34| ). Taking the limit s — > em- 
ploying ( |3.10| ) we obtain the TQ-equation, 

Q+(z)t (1 \z) = q a - SA T^(zq)Q + (zq- 2 ) + q S ^ a T^ {zq- l )Q+ {zq 2 ) . (3.17) 

A similar relation holds for Q~ when employing the transformations ( |2.27| ), ( p.30| ) and 
(3.12). The TQ equation holds true also for vanishing external magnetic field, i.e. a = 0. 



If the magnetic field is nonzero, however, there is another identity which makes use of both 
solutions and on which we will focus. It yields a simple expression for all elements in 
the fusion hierarchy in terms of Q^, 

T (d-l) (z) _ q d{a ~ SA) Q + (zq- d )Q-(zq d ) - q d ( S ^Q + (zq d )Q-(zq- d ) 

\ > qa-SA — qS A -a ' \ ■ J 

The last expression can even be analytically continued to complex dimension d if one wishes 
to make contact with the trace functional introduced in the context of correlation functions 
for the infinite XXZ spin-chain. Here we shall not pursue this aspect further but refer the 



reader to [11] and references therein; see also 



4. The Wronskian relation 



Note that when setting d = 1 in ( 3.18| ) the left hand side of the above equation is explicitly 
known and we arrive at 

- m,fa - ir = -_ ywg-cr') . (4 .i) 

This functional relation, which is the discrete analogue of a Wronskian in the theory of 
second order differential equations, can therefore be employed to compute the eigenvalues 
of Q^. It is this observation which we will investigate further in the remainder of this 
paper. 
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4.1 System of quadratic equations 



Expanding the Wronskian relation (4.1) with respect to the spectral parameter z yields 
the following system of iV-quadratic equations for the unknown coefficients e. which are 
the elementary symmetric polynomials in the Bethe roots, 

V- ( n \ ( n \ < sk-l V- q a - SA ~ k+l - q SA-a+k-l 

£ UXJ M ^ q «s A - q s A - a 't*i ■ (4-2) 

k+l=m V / \ / k+l=m H * 

It should be emphasized once more that this set of equations contains for a ^ all the nec- 
essary information about the spectrum of the quantum transfer matrix as well as the higher 
spin transfer matrices. Recall that N is the Trotter number, w = e~@ ^ N = exp(— — -) 
contains the temperature variable and q a = e h ^l 2 the external magnetic field. Setting now 
n = 2 in ( |3.18| ) the eigenvalues of the quantum transfer matrix Ql.2| ) are obtained from the 
identity 

, , _ q 2( - a - S ^Q + (zq- 2 )Q-(zq 2 ) - q 2 ^-*)Q+( zq 2)Q-(zq- 2 ) 
T ^ Z,W) {q a ~ s A - q S A -a^ zwq _ q -iY( z /wq- q) n ' 

The equations ( |4.2| ) and ( ^4.3| ) are the aforementioned generalisations of the identities ( 1.15| ) 
and ( |l,16| ) in the introduction. 

4.1.1 The Bethe ansatz equations 

Having identified ( [4,2; ) as the key relation in describing the spectrum of the quantum 
transfer matrix for finite N we need to discuss the relation with the Bethe ansatz equations 
which are usually considered to be the fundamental set of identities for discussing the 
spectrum. Starting from the quantum Wronskian we set 



/+ -q SA - a Q+(q 2 /xt)Q-(l/xf) 2 + + 



Z = 1 l X i- q a-S A _ q S A -a = ( W / X i ~ 1) (<? / WX i ~ 1) 

and obtain 

_ q -2« fr ±M_ ~ I' 1 = ( W q-xtq^ \ n ( l-wxf \ n (44) 

The last equation is one among the n + Bethe ansatz equations as they can be found for 
instance in ||]; see equations (29) and (30). To facilitate the comparison note that under 
the parametrisation 

x i = e lX \ q = e i ' f , w = e~^ T = e- {q ~ q ' 1)f3/N , q 2a = e Ph 

the Bethe ansatz equations are rewritten as 



<j)(\j + i) sinh 2 l(Aj - ir + 2%) sinh 2 |(Aj + ir) 
<K A j ~ i) ' sinhf + ir - 2») sinhf ^(Aj - ir) 



o( 3h sinh^Aj - X k + 2i) _ <>jSft Q(Aj + 2i) 



A sinh ^(Aj - A fc - 2i) ' " C Q(A J - - 2«) 
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which is the notation used in |2j. From this we infer that the Wronskian relation ( |4.1| ) 
implies the Bethe ansatz equations, the converse is not true. Note that Q~ yields another 
set of Bethe roots, which in terms of the algebraic Bethe ansatz (see the appendix) allow 
one to construct the eigenvector from the lowest (instead of the highest) weight vector. 

4.2 Special solutions in the Sa = sector 

As pointed out earlier the quantum transfer matrix block decomposes with respect to the 
alternating spin-operator Sa- In particular, the largest eigenvalue Ajy appears to be in the 
Sa = sector. Since limjv^oo An yields the partition function (|l.3|) in the thermodynamic 



limit its corresponding solutions to (4J2) are of particular, physical interest. 

Quite generally, we infer from ( [4.1|) that there are certain symmetries among the so- 
lutions to the Wronskian relation, for instance when replacing z — > z . It is then natural 
to assume that some solutions are invariant under these transformations, especially when 
they belong to non-degenerate or distinguished eigenvalues such as Ajy. Thus, one would 
expect that there is a subset of solutions for which 

Q ± (z)=f[(l-z/xf) (4.5) 
i=l 

holds. That is, are the inverse or reciprocal polynomial of each other. In fact, one 
verifies numerically for Trotter numbers up to 16 that in the Sa = sector there are 
always 2 n such "special" solutions (for q real or on the unit circle) and that among this 
set of solutions is the one describing the largest eigenvalue of the quantum transfer matrix. 
Conjecturing this to be true for all N one can then halve the number of variables in ( |1.2| ) 
arriving at the result ( 1.15| ), ( |1.16| ) stated in the introduction. A similar observation has 



been made previously for the twisted XXX spin-chain; see Section 5.1 in [15]. 

Provided that q (and therefore w = exp[(g — q~ 1 )(3/N]) lies on the unit circle there is 
another obvious "symmetry" of ([4.1|), complex conjugation. In fact, one finds for Sa = 
that there exist solutions invariant under this transformation, i.e. obeying the addi- 
tional restriction 

n 

Q ± (z) = H(l-(xfyz), zeR. (4.6) 
i=i 

In terms of the elementary symmetric polynomials the restrictions (|4.5| ) and ( [4.6[) cor- 
respond to the identities 

e* and 4 = (et)*> ( 4 - 7 ) 

respectively. Whether these solutions persist for Trotter numbers N > 16 needs to be 
investigated numerically. We leave this to future work. 

5. Conclusions and outlook 

In this paper Baxter's Q-operator has been constructed for the quantum transfer matrix 
of the XXZ spin-chain. The main motivation has been to derive a set of equations which 
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allow one to describe the spectrum of the quantum transfer matrix and are simpler than 
the previously known Bethe ansatz equations. For non- vanishing external magnetic field 
this is indeed possible and the order of the equations can be reduced from N (Bethe ansatz 
equations (4.4)) to order two (the Wronskian relation Q4.2| )) which is a great simplification 
for numerical computations. Of particular interest is the largest eigenvalue of the quantum 
transfer matrix which in the thermodynamic limit (i.e. the physical system size L tends to 
infinity) contains all the relevant finite temperature behaviour of the spin-chain. For this 
eigenvalue we observed further simplifications in its polynomial structure which lead to a 
reduction of the number of variables by a factor two. The latter observation is based on 
symmetries of the Wronskian relation ([O]) and numerical computations which were carried 
out up to the Trotter number N = 16. To go beyond this bound requires more extensive 
numerical computations which are planned to be carried out in future work. 

Another aspect which has been omitted from the present work is the derivation of 
integral equations similar to those non-linear integral equations which have been previously 
obtained on the basis of the Bethe ansatz equations; see e.g. P] and references therein. 
The basis for such a derivation is numerical evidence for the distribution of the Bethe 
roots in the large Trotter number limit, N — > oo. For this reason the discussion of integral 
equations and their connection with the thermodynamic Bethe ansatz are postponed until 
the necessary numerical data has been obtained. 

For vanishing magnetic field the Wronskian relation ceases to be valid but in the text 
it was shown that at roots of unity one might be able to employ the same loop algebra 
symmetry as the one which exists for the six-vertex model. This might possibly help to 
reduce the order of equations or simplify the computation of the spectrum of the quantum 
transfer matrix. 
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A. Properties of Q for generic q 

For comparison with the line of argument for the classical six-vertex model the reader 
might wish to consult ||. The line of argument closely follows the exposition given there. 

A.l The fusion hierarchy in terms of Q 

Setting the free parameter s in p + to the special value s = q 2d , d = 1,2,... one derives 
the following identities by restricting the £ and £*-operator given in Q3.4| ) and ( |3.5| ) to the 
subspaces V <( i = s P an {|^)}fc=o an d V>d = span{|/c)}^ rf , 

s = q 2d : £?{z)\ v<d = q d (l <g> ^L*^ 1 ) ( zq d )(l <g> q-***) (A.l) 
s = q 2d : £*(z)\ v > d = q 2d (l ® q dr7Z )£*(zq 2d ; s q- 2d )(l ® q~ 2duZ ) (A.2) 

Employing that by construction the Q-operator commutes with the alternating spin-operator, 
[Q, Sa] = 0, one deduces from these last two equations the following identity for the higher 
spin quantum transfer matrices, 

T^iz) = q di ~ a ~ S ^Q{zq~ d - s = q 2d ) - q d{SA - a) Q(zq d ; s = q~ 2d ) . (A.3) 

Below we will see that this relation can be simplified further due to a factorisation of the 
Q-operator into two parts. 

A.2 Algebraic Bethe ansatz computation 

In order to obtain the eigenvalues we apply the algebraic Bethe ansatz for the quantum 
transfer matrix (see e.g. 0) and compute the action of Q on a Bethe state. Denote by 
{1)1} by the orthogonal basis in C 2 and consider the reference state 

|o> =1 ® T • • • ® J. ® T • ( A - 4 ) 

Let {A, B, C, D} be the matrix elements of the quantum monodromy matrix ( |2.21| ) then a 
Bethe state is given by 

x+,...,x+ + ) :=B(l/x+)---B(l/x+ + )\0) 

with {x+ , —,x% } being a solution to the Bethe ansatz equations ( f4.4[ ). Since the auxil- 
iary space of the quantum transfer matrix (|L^) and Q-operator ([T?]) are the same as in 
the ordinary case of the classical six-vertex transfer matrix with quasi-periodic boundary 
conditions the results from Q apply. Using the commutation relations between the matrix 
elements A, B, C, D and those of the monodromy matrix of the Q-operator, 

Q kjl (z; s) := (k\ q ah ^^ N {zw)Sl* N ^{z/w) ■ ■ ■ £ 2 (zw)£l(z/w) \l) , (A.5) 

detailed in H we obtain 

Q(z; s) = 

2+ (A: + 1| Oj \k + 1} (k\ dj \k) - (k\ Cjbj \k) " 

.k=0 



x l ) •••) x n + 
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Here we have introduce the abbreviations 

(U) (A - 6) 

for the matrix elements of the £-operator. Inserting the explicit expressions for the latter 



which can be read off from (3.4) one arrives at 
Q(z; s) \x~^ , x^ 



-2k-l 



(l-zx+)(l 



^2 (0| Q kk (z; s) |0) n ^"+777 

j=i (1-zsq 2k xj)(l 



ZS XI 



zsq~ 



-2fc-2 



X 



1 ' 



n + 



q S *- a Q + ( Z )Q + ( Z s)J2 



7 2fc(S A -a)_(0) 



(zsg" 



fc=0 



Q+(zsg- 2fc )Q+(zsg- 2fc - 2 ) 



In the last line we have used the definition of Q + as a polynomial Q3.14| ) and 

(0| Q kk (z;s) |0> = q^ 2k+1 \(zwsq- 2k - l)(zsq- 2k - 2 /w - l)f • 

Taking the limit s — > in order to fix the normalisation constant we deduce from this 
expression the following formula for Q~, 



Q-(z) = Q- 



fc=0 



2k(S A -a) T (0)f zq -2k-l\ 
Q+(zq- 2k )Q+(zq- 2k - 2 ) ' 



(A.7) 



Note that in ||] the vanishing of the unwanted terms in the Bethe ansatz has only been 
verified for n+ = 1,2,3, since the algebraic Bethe ansatz computation is more involved 
than in the case of the transfer matrix due to the infinite-dimensional auxiliary space of 
the Q-operator. However, the result for the spectrum coincides with the one found at 
roots of unity and the eigenvalues satisfy all aforementioned functional relations for the 
Q-operator, which have been derived by different means. We shall take this as sufficient 
evidence that the algebraic Bethe ansatz computation presented above holds also true for 
n 4- > 3. 



B. Properties of Q when q is a root of unity 

Let us now turn to the case when q is a primitive root of unity of order I. 
B.l Functional relations for the fusion hierarchy 

The evaluation module p + is now finite-dimensional according to Q3.3| ). For analysing 
the spectrum of Q we now rely on a functional relation derived in B. To connect with 



the discussion therein we set r = and s = /i in ( |3.l| ) and obtain the evaluation 
representation iTz specified in equation (15-16) of M. Thus, in the following we refer to 
/i~ 2 ) as iTz- Then the following short exact sequence holds (see equations 



p^[r = fi ,s 
(52-53) in i) ; 



7T 



jivq 

m 



A 1 

7T 2 



7Ti 



7T 



fivq 

m 



-e'+i 



<8> 7T 



vq 



e'+i 



. 



(B.l) 
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(£' —2) 

Here tt z is the evaluation representation of spin {£ — 1)/2. As before this decomposition 
of the tensor product of representations implies a functional relation, but this time it 
involves a product of two Q-operators with different values for the free parameters entering 
Q3.1|). Setting s = pT 2 and t = v~ 2 this functional relation reads 



Q(zq 2 /s;s)Q(z;t) = q SA ~ a Q{zq 2 / s; stq- 2 )[{zq 2 - l) n (z - If + /(S A -a) T (e'-2) , 

(B.2) 

compare with equation (46) in |j[c|| . Assume that [Q(z\; s\), Q(z2', S2)] = for arbitrary 
pairs zi 5 2,si,2 G C. This has been proved for £ = 3,4,6 by explicitly constructing the 



corresponding intertwiner for p + (z\,si) <g) p (#2; ^2) [|1G|, |17|| . Thus, the eigenvalues of 



Q(^; s) must be polynomial in z (and s) their most general form being 

1 _ q 2£'{S A -a) k N-k 
<2(^ S ) = n a-S A _ n S A -a I! (1 " Z X 3 ) I] (1 " Z yj(s)) . 

h q 3=1 j=i 

Here we have taken the limit z — > to fix the normalisation constant, 

lim Q(z; s) = Tr q ( a - s ^i _ ( B-3 ) 



2^0 



We assume that the roots Xj = x~j are independent of s while the yj's depend on it allowing 
for the possibilities that either k = or k = N. Inserting this general expression for the 



eigenvalue into (B.2) one deduces that the roots yj(s) can only depend linearly on s, i.e. 
Uj(s) = xjs for some xj . This implies the factorisation ( |3.9| ) in the text. Futhermore we 
infer from (3.12) that k = N/2 — Sa in each fixed Sa sector. 
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